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ABSTRACT 

In the last ten years, observational evidences about a possible connection between 
Type la Supernovae (SNe) properties and the environment where they explode have 
been steadily growing. In this paper I discuss, from a theoretical point of view but with 
an observer's perspective, the usage of hght echoes (LEs) to probe the CSM around 
SNe of Type la since, in principle, they give us a unique opportunity of getting a three- 
dimensional description of the SN environment. In turn, this can be used to check the 
often suggested association of some la's with dusty/star forming regions, which would 
point to a young population for the progenitors. After giving a brief introduction to the 
LE phenomenon in single scattering approximation, I derive analytical and numerical 
solutions for the optical light and colour curves for a few simple dust geometries. A 
fully 3D multiple scattering treatment has also been implemented in a Monte Carlo 
code, which I have used to investigate the effects of multiple scattering. In particular, 
I have explored in detail the LE colour dependency from time and dust distribution, 
since this is a promising tool to determine the dust density and derive the effective 
presence of multiple scattering from the observed properties. Finally, again by means 
of Monte Carlo simulations, I have studied the effects of multiple scattering on the 
LE linear polarization, analyzing the dependencies from the dust parameters and 
geometry. Both the analytical formalism and MC codes described in this paper can 
be used for any LE for which the light curve of the central source is known. 

Key words: supernovae: general - ISM: dust, extinction - ISM: reflection nebulae - 
methods: analytical - methods: numerical - radiative transfer - polarization 



1 INTRODUCTION 

In the last few years the study of light echoes (hereafter 
LEs) in Supernovae (SNe) has become rather fashionable, 
since they provide a potential tool to perform a detailed to- 
mography of the SN environment (see |&iEorman 2003 and 
references therein) and, in turn, can give important insights 
into the progenitor's nature, a matter which is still under 
debate. Due to the typical number density of dust parti- 
cles which are responsible for the light scattering, LEs are 
expected to have an integrated brightness about ten magni- 
tudes fainter than the SN at maximum. For this reason, 
a SN la in the Virgo cluster is supposed to produce, if 
any, an echo at a magnitude V~21.0. This has the sim- 
ple consequence that it is much easier to observe such a 
phenomenon in a la than in any other SN type, due to 
its high intrinsic luminosity. As a mat ter of fact, only four 
cases are known: t he SNe la 199 1T JSchmidt et al.|ll994l 
ISnarks et all 119991) and 1998bu iCappellaro et alJ l200ll). 
and the type II SN1987A jCrotts. Kunkel fc McCarthy! 



19891 ICrotts fc Kunkel Il99ll . IXu Crotts fc KunkeJ Il994l 
Crotts. Kunkel fc Heathcotd Il995l. IXu. Crotts fc KunkeJ 
19951) and 1993 J JSueerman fc Crottjl2002t) . As expected, 
the LE detections for the two core-collapse events occurred 
in nearby galaxies: LMC (d=50 Kpc) and M81 (d=3.6 Mpc) 
respectively. Now, while a dusty environment is expected 
for explosions arising in massive and short-lived progenitors, 
this is in principle not the case for the kind of stars which 
are commonly supposed to generate la events, i.e. old popu- 
lation and low mass. This is why the detection of substantial 
dust in the immediate surroundings of thermonuclear SNe, 
or at least in a fraction of them, would indicate a different 
scenario. 

Actually, several authors have pointed out that the 
observed characteristics of la's, such as intrinsic lumi- 
nosity, colour, decline rate, expansion velocity and so 
on, appear to be related to the morphological type of 
the h o st galaxy jFil^^enkol^si jBraiwl^^^ai^BnBergD 
| l993l Ivan den Bergh fc PazdeJ 19921: iHamuv et al.l Il99a 
iHamuv et alJl2000l: lHowellll200ll) . Since these objects rep- 
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resent a fundame ntal tool in Cosmology (see for instance 
lLeibundgutl2000h . it is clear that a full understanding of the 
underlying physics is mandatory in order to exclude possible 
biases in the determination of cosmological parameters like 
n and A. In this framework, disentangling between possible 
la sub-classes is a fundamental step. In this respect, an im- 
portant year in the SN history is 19 91, when two extreme ob - 
jects were disco vered, i.e. SN1991T ( Filipp enko et alll992al) 
and SN1991bg iFilippenko et al. 1992b'l. The former was an 
intrinsically blue, slow declining and spectroscopically pecu- 
liar event, while the latter was intrinsically red, fast declin- 
ing and also showing some spectral peculiarities. From that 
time on, several other objects sharing the characteristics of 
one or the other event were discovered, indicating that these 
deviations from the standard la were, after all, not so rare. 
Of course, one of the most important issues which was gener- 
ated by the discovery of such theme variations concerned the 
explosion mechanism and, in turn, the progenitor's nature. 

The growing evidences produced by the observations 
in the last ten years have clearly demonstrated that the 
sub-luminous events (1991bg-like) are preferentially found 
in early type galaxies (E/SO), while the super-luminous ones 
(1991T-like) tend to occur in spirals (Sbc or later). This 
has an immediate consequence on the progenitors, in the 
sense that sub-luminous events appear to arise from an old 
population while super- luminous ones would rather occur 
in star-forming environments and therefore would be asso- 
ciated with a younger population. This important topic has 
been recently discussed in a work bv lHowelll J200j) . to which 
I refer the reader for a more detailed review. What is im- 
portant to emphasize here is that 1991T-like events tend to 
be associated with young environments and are, therefore, 
the most promising candidates for the study of LEs. Or, 
in turn, if LEs are detected around such kind of SNe, this 
would strengthen their association with sites of relatively 
recent star formation. 

The first case of LE in a la (1991T) seems to confirm 
this hypothesis, in the sense that the SN was over-luminous 
and the host galaxy (NGC4527) is an Sbc and also a liner. 
Slightly less convincing is the other known case (1998bu), 
since the galaxy (NGC3368) is both an Sbc and a liner, but 
the SN is not spectroscopically peculiar. The only charac- 
teristic in common with SN1991T is its decline rate Amis, 
which is lower than average, even though not so extreme 
as in th e case of 1991T. Never theless, the HST observa- 
tions by iGarnavich et alJ (|20o3)) show that a significant 
amount of dust must be present within lOpc from the SN. 
Of course no statistically significant conclusion can be drawn 
from such a small sample, which definitely needs to be en- 
larged. For this reason, during the past years, I have been 
looking for new cases, the most promising of which was rep- 
resented by SN1998es in N GC 632. This SN , in fact, was 
classified as a 1991T-like bv ljha et all (Il998l) . who also no- 
ticed that the parent galaxy was an SO, hosting a nuclear 
starburst JPogge fc Eskridedliggl . Moreover, the SN was 
found to be projected very close to a star forming region 
and to be affected by a strong reddening, which all together 
made SN 1998es a very good candidate for a LE study. 

While the present paper is devoted to the theoretical as- 
pects of the LE phenomenon, in a forthcoming article (Patat 
et al. in preparation) we will be dealing with the application 
of these results to the known cases of SNe 1991T, 1998bu 



and the test case of SN 1998es. In the same work we will also 
present high resolution spectroscopy for the three events and 
other unpublished data relevant for the LE discussion. 

The paper is organized as follows. After giving a gen- 
eral introduction to the LE phenomenon in single scattering 
approximation in Sec. |21 in Sec.|3l derive analytical solu- 
tions for a thin perpendicular sheet and a spherical shell 
to illustrate the effects of forward scattering. In Sec. 0| I 
then present a numerical solution for a double exponential 
dust distribution, typical of a spiral galaxy. Single scattering 
approximation is refined in Sec. 13 where I introduce the at- 
tenuation and the concept of LE effective optical depth. The 
inclusion of multiple scattering is presented in Sec.|Sl where 
I illustrate the implementation of a Monte Carlo (MC) code, 
whose results are shown and discussed in Sec. Q The same 
code is then used to compute the LE spectrum (SeclSJ and 
broad band (B — V) colour (Sec.|^, which is compared with 
the results obtained with the single scattering plus atten- 
uation approximation. The effects of a change in the dust 
mixture are described and discussed in Sec. 1101 In Sec. 1111 
I introduce the concepts of MC polarization calculations, 
which I have used in a code to study the effects of multiple 
scattering, presented in the same Section. Finally, in Sec. 1121 
I give a closing discussion and in Sec. ll3I I briefly summarize 
my results. 



2 SINGLE SCATTERING APPROXIMATION 

Starting with the pioneering work by ICouderj il939l) . 
the problem has been addressed by several au- 

thors (s ee fo r example |p wek 1983, Chevalier 198^, 

Schaefeilll987l. lEmmering fc Chevalier 1989, Spark^^M 
Xu. Crotts fc Kunkell Il994l . iSparks, .1997l and iSugermanI 
2003h . who have all adopted the single s cattering approx i- 
mation. The only exception is the work bv lChevalieil il983l . 
who has used a MC code to study the case of a spherically 
symmetric stellar wind with a dust density law. 

Given the number of existing publications on this topic, 
I will give here only a brief introduction to the LE phe- 
nomenon in the single scattering approximation, while I will 
focus on the analytical and numerical solutions which I shall 
use later on to test the MC code. 

2.1 Integrated echo light curve 

Let us imagine a SN immersed in a dusty medium at a dis- 
tance d from the observer. If we then consider the SN event 
as a radiation flash, whose duration AtsN is so small that 
c AtsN <^ d, at any given time the LE generated by the 
SN light scattered into the line of sight and seen by the ob- 
server is confined in an ellipsoid, whose foci coincide with 
the observer and the SN itself. Since d is supposed to be very 
large compared to all other geometrical distances, this ellip- 
soid can be very well approximated by a paraboloid, wit h 
the SN lying on its focus (see for example IChevali'eilll986ll . 
More generally, the paraboloid can be regarded as the locus 
of those scattering elements which produce a constant delay 
t and for this reason I will refer to it as the iso-delay surface. 
If L{t) is the number of photons emitted per unit time by 
the source at a given wavelength, the flux of scattered pho- 
tons per unit time and unit area which reach the observer 
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Figure 1. Geometry of the problem in single scattering approx- 
imation. 



at time t is obtained integrating the flux scattered by the 
volume element dV (see Fig.0 over all iso-delay surfaces: 



47rd^ F{t) = / L{t- t') f{t') dt' 



(1) 



The f{t) function, which has the dimension of the in- 
verse of a time, contains all the information about physical 
and geometrical properties of the dust and it is defined as 
follows: 



f{t) = Cext UJ C 



m 



ct/2 



n(tf, z, i) dip dz 



(2) 



where ip is defined as in Fig. ^ r — z + ct is the dis- 
tance of the volume element from the SN, n = n{x, y, z) 
is the number density of scattering particles, Cext is the 
extinction cross-section^, lo is the dust albedo, 9 is the 
scattering angle defined by cos 6 = z/r and ^(6) is the 
scattering phase function, normalized in order to have 
^,_^(S)dQ. = 1. As usual in this kind of studies (cfr. 
[chevalier 198^ . I have adopted the formulation proposed 
by Henyey & Greenstein (1941, hereafter H-G), which in- 
cludes the degree of forward scattering through the parame- 
ter g: isotropic scattering is obtained for g=Q ($(0) = l/47r) 
while g=l gives a complete forward scattering. Both em- 
pirical estimates and numerical calculations for the optical 
wavelength range give g ~0.6 l|Whitelll979l) . which corre- 
sponds to an average scattering angle 9 ~50°. More sophis- 
ticated calculations, including a full Mie's treatment, show 
that the H-G formulation tends to underestimate the for- 
ward scattering, but the approximation is reasonably good 
(see for example iBianchi. Ferrara fc Giovanardil)l996l) . As 



for the dust albedo, this is usual ly assumed to be u> ~0.6 
jMathis. Rumpl fc Nordsieldll977t) . 

From Eq. Q it is clear that the out-coming scattered 
signal is the convolution of the SN signal L{t) with the im- 
pulse response function fit) and therefore, at any given time 
t, the observer will receive a combined signal, which contains 
a sum of photons emitted by the SN in the whole time range 
— t. This equation is formally identical to the one which de- 
scribes the resulting acoustic wave reverb eration in a given 
environment (see for example ISpiutiEoOlll . 

The global flux received by the observer is obtained 
adding to the echo contribution the photons coming directly 
from the SN, after correcting for possible extinction: 



Lit) 



47rd2 
where 



+ Fit) 



Td 



Cext n(0, 0, z) dz 



(3) 



(4) 



is the dust optical depth along the line of sight. Of 
course, if other dust is present but it does not contribute 
to the LE, because it is too far from the source, both the 
SN and the LE would suffer by additional extinction. 

Due to the very short duration of the SN burst and 
in the single scattering approximation, there is a very good 
correspondence between the observed impact parameter p = 
ct i2z + ct) and the distance r of the scattering material 
from the SN (see Fig.[TJ. This can be easily calculated as 



ct 



+ ct 



(5) 



and it offers the unique possibility of measuring the 
third dimension of a resolved LE image with a spatial reso- 
lution of the order of cAtsN, provided that d is known with 
some approximation. I note that this is always true, no mat- 
ter what the dust geometry is, only provided that the dust 
optical depth is not too large to cause multiple scattering, 
which would introduce some degree of degeneracy in the p 
vs. r relation (see Sec.0. 



2.2 Echo surface brightness 

If, for a given point (a;, y) on the sky plane, we consider 
the sum of all contributions along the z axis, we can easily 
compute the echo surface brightness E at any given time as: 



Eix,y,t) = 



2c LO Ce. 

47rd2 



L(t-t') 



nix,y,t) $ie) dt' (6) 



The explicit dependency of $ from x, y and t through 
the scattering angle 9 can be obtained remembering that 
cos 9 = z/iz + c t) which, using the paraboloid equation (cfr. 
IChevalienll986ll . can be reformulated as: 



cos9ip,t)=[ 



cH'^ 



+ 1 



(7) 



^ It must be noticed that in the literature Cext is usually given 
per Hydrogen atom (see for example Drainc 2003). In that case, 
n has to be replaced with the Hydrogen numerical density n(i?). 



3 SIMPLE ANALYTICAL SOLUTIONS 

Eqs.|lland[Tl can be solved numerically, using the observed 
light curve of a typical la for Lit). Nevertheless, instructive 
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results can be obtained analytically if one assumes that the 
SN light curve is a flash: 



Lit) 



Lq for t ^ AtsN 
for t>AtsN 



(8) 



In this approximation and for t ^ AtsN, the convolu- 
tion in Eq.0 becomes trivial and gives: 



47V(t Fit) = Lo AtsN fit) 



(9) 



and hence the whole problem reduces to compute fit) 
defined by Eq. El The flash durat ion can be estimated from 
the observed light curve as Lo AtsN = °° Lit) dt, and 
turns out to be of the order of 0.1 years for a typical la. The 
use of H-G formulation with g ^0 in Eq. |21 makes the an- 
alytical integ ration possible o nly under certain conditions. 
For example. IChevalieil il986l) has found an analytical solu- 
tion for the case of a density law n oc r~ , typical of stellar 
winds. Here I will focus on two simpler cases only, those of 
a thin perpendicular sheet and a thin spherical shell, which 
are useful to illustrate the effect of forward scattering. In 
both geometries the solution becomes very easy, since p can 
be considered constant along the iso-delay surfaces where 
n /O. 

Let me first consider a thin sheet perpendicular to the 
line of sight, placed at a distance R from the SN and with 
a thickness AR. A limit on AT? in order to fulfill the above 
condition can be derived differentiating the relation between 
H and z (see Sec.|2J and imposing <^ 1. This gives: 



Ai? < 



iR + ct)^ 
ct 



Having this in mind and recalling that for z = R, one 
has that /i — z/r = R/iR + ct), Eq.|5|can be integrated and 
it gives: 



fit) 



1^ Cext c no ( 1 - g"^) 
2(1+52-25 



^i^lnfl + 



AR 

R + ct 



(10) 



where, with the additional assumption that AT? <ti R + 
ct (not necessarily included in the previous condition), the 
logarithmic term can be app r oxima ted with AT?/(T? -I- ct), 
as done bv ICappellaro etaO i200J) . As for the integrated 
luminosity, I notice that under these circumstances the time 
dependency is very mild. Therefore, using typical values for 
the relevant parameters (5=0.6, uj—0.6, Atsiv=0.1 yrs) and 
substituting Eg. llUl in Eq. |31 the normalized LE flux can be 
written as 



Fit) 



0.3 



Td 



Lo R 

where I have posed = Cext no AR (see Eq. 2J ■ This 
basically means that, at least for a distant cloud (T? 2> ct), 
the LE luminosity is inversely proportional to R and di- 
rectly proportional to r^. For example, for a sheet placed 
at T?=100 lyrs from the SN, in order to produce a LE 10 
mags fainter than the SN at maximum, one needs an optical 
depth Td ~0.03. 

In general , the resulti ng echo will appear as a ring of 
radius pr = rf(2T? -I- ct) and thickness ApR ~ ct AR/ pR, 
with roughly constant surface brightness. 

The overall effect of forward scattering can be evaluated 
directly from Eq.ITUl In fact, for R > ct, it is R/iR + ct) ~1, 
so that the forward scattering term becomes approximately 



(1 + 5)/(l — 5)^. This means that going from 5 = to 
5 = 0.6, the echo maximum luminosity increases by more 
than a factor of 10. Moreover, as a consequence of the loss 
in the scattering efficiency due to the increase of the scatter- 
ing angle with time, the resulting light curves get steeper. 
This is particularly true when the dust is not too far from 
the central source, while for increasing values of R the time 
dependency of the forward scattering term tends to disap- 
pear (see Eg. lion . 

The inclusion of forward scattering produces also a sen- 
sible flattening of the surface brightness profile. In fact, an 
increase of 5 causes the material which is close to the SN to 
be less efficient in scattering the incoming photons, and this 
tends to balance with the fact that more distant material 
receives less photons from the SN but it is more effective in 
deflecting its light. Furthermore, for homogeneous distribu- 
tions, the surface brightness at p = ct decreases by a factor 
(1 — 5'^)/(l -1-5'^)^''^ which, for 5 = 0.6, corresponds to about 
one magnitude. This is relev ant for the de tection of the max- 
imum polarization ring (see lSparkj|l994l and the discussion 
in Sec. Illl here'l. 

Another example of analytical solution in the case of 
5 7^0 is the one of a thin spherical shell of radius R and 
thickness AT? <C T?. In this case, since R = z + ct, we have 
that p — iR ~ ct) /R and therefore the solution of Eq. |5| 
gives: 



fit) 



Cext c no (1 - 5^) 
2[l+g2-2g^^^fr2 



In 



Also in this case the LE would appear as a ring, whose 
radius pr = ct(2T? — ct) grows, reaches a maximum at 
t = T?/c, decreases again and finally shrinks to zero for 
t — 2T?/c. As in the case of the thin sheet, the echo lumi- 
nosity is approximately proportional to Td R~^ and forward 
scattering has a similar effect. 



4 NUMERICAL SOLUTION FOR A DOUBLE 
EXPONENTIAL DUST DISK 

The numerical integration of Eq. |5| allows one to solve 
the problem under the single scattering approximation for 
whichever dust geometry and input SN light curve. 

For this purpose I have chosen to adopt a V tem- 
plate light curve obta ined smoothing t he comb ined data 



sets of SNe 1994D Patat et al.) Il99d iRichmond ct alJ 

E ICappellaro et ali Il99'i1) and 1992A ISuntzeff .199^ 
idlar^^^n^g^ . t wo very similar an d well sampled 
objects (see for example IPatat et al.lll99^ . The result is 
shown in Fig. |5| where I have plotted the corresponding 
(T? — V) colour curve as well. The latter will serve to com- 
pare the intrinsic SN colour with that of the LE (see Secl^J. 

In Fig. 121 I have traced the results obtained for a 
exponential face-on dust disk, whose density profile was 
modeled using the typical expres sion (see for example 
iBianchi. Ferrara &i Giovanardilll996h : 



i(T?, Z)=no ^^p(~"^) ™P 



Zd 



where T?, Z are the cylindric galactic coordinates and 
T?£j, Zd are the characteristic radial and vertical scales. Fol- 
lowing Bianchi et al. (1996), I have adopted T?d=4.0 kpc and 
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Time since explosion (yeors) 



Figure 2. Template light curves for B (black) and V (gray) 
passbands. The solid lines trace the smoothed version used in the 
calculations while the filled dots indicate the observed data. The 
lower panel shows the corresponding {B — V) colour curve. 

Zd~0.14: kpc, which can be considered as typical values for 
a spiral galaxy. The density no is constrained by the cen- 
tral optical depth of the disk seen face-on, which is given by 
t(0) = 2nn C^xt Zd - Imposing a typical value t{0) — 1 (see 
IXilouris etall 1 1999ft and for Cext~ 5xl0'^'^ cm^, no turns 
out to be 2.3 cm"^. The SN was placed at i?sjv=2xl0* lyrs 
(~6 kpc) and Zgjv=300 lyr above the galactic plane, which 
imply an optical depth Td ~0.06. 

For both isotropic and forward scattering I have run the 
calculations using a flash (dashed line) and the template V 
light curve (thick solid line) and this shows that the results 
are pretty similar, except for the early phases, where the 
flash tends to produce a sharper and brighter peak. This is 
anyway irrelevant, since at such epochs the SN is typically 
more than 10 magnitudes brighter that the LE and therefore, 
at least for la's, Eq.|5]gives a fair description of the real light 
curve. 

The effect of forward scattering is clearly illustrated in 
Fig. 13 with g=Q.& the light curve starts to deviate from the 
radioactive tail when the SN is about 1.5 mag brighter than 
with g=Q. Moreover, forward scattering has an effect on the 
resulting slope for t >3 yrs, which is 0.1 and 0.05 mag yr~^ 
in the two cases respectively. 

In conclusion, under rather normal conditions, a Type 
la SN exploding in the disk of a spiral should always produce 
a LE which is of the order of 10 mag fainter than the SN 
at maximum, without the SN being heavily reddened. Of 
course, placing the SN in the inner parts of the disk or on 
the far side of the galaxy would increase the echo luminosity 
and the extinction suffered by the SN itself. For example, 
leaving all the other parameters unchanged and placing the 
SN at ZsN=0 would enhance the LE by 0.7 mag, while the 
optical depth would grow to Td ~0.1, which is anyway still 
a rather low value. Since there are good reasons to believe 
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Figure 3. Single scattering numerical solutions for a face-on ex- 
ponential disk for 3=0 (lower curves) and 9=0.6 (upper curves). 
The dotted curves trace the global light curves (SN-|-echo) while 
the long dashed line indicates the radioactive tail of the tem- 
plate light curve. In all examples Cext=5xl0~^^ cm^, aj=0.6, 
_R5]V=2xlO'* lyr and Zsn=300 lyr were used (see text for more 
details). 

that 1 ^ r(0) ^ 5 iBianchi |2003) . a Type la within 1-2 
dust scale heights should always produce an observable LE, 
unless it is located very far from the galaxy center. 

Using numerical integration, I have explored a few other 
simple cases, like front sheets with different inclinations, thin 
spherical shells, winds and spherical clouds placed at dif- 
ferent distances from the SN and with varying offsets with 
respect to the line of sight. In the latter case, when there is 
no scattering material along the line of sight, the LE starts 
to take place at t >0 and, in this respect, the parallel with 
acoustic physics becomes more pertinent, in the sense that 
the reflected signal starts to be detached from the direct 
impulse. Actually, in all cases where there is material on 
the line of sight, one should rather talk about light reverbs 
(LRs), which corresponds more appropriately to the analo- 
gous effect for acoustic waves. 

In general, it turns out that all dust geometries produce 
similar LE luminosities and light curves. For typical values 
of no, Cext, and g, the LE at three years after the explosion 
has an integrated brightness between 9 and 11 mag below 
SN maximum. Therefore, unless one is able to resolve the 
LE image, it is very difficult to disentangle between different 
dust geometries on the basis of the light curve shape alone. 



5 SINGLE SCATTERING PLUS 
ATTENUATION 

The single scattering approximation can be refined in order 
to take into account the presence of absorption. Some times 
this is referred to as single scattering plus attenuation (here- 
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to the observer 

A 



Table 1. Relevant parameters for the test-case geometries. 




Figure 4. Meaning of the optical path I in the single scattering 
plus attenuation approximation. 



after SSA. See for example I Wood et alJIlQQ^) . This approx- 
imation is obtained assuming that the fraction of a photon 
packet which reaches the observer has undergone one single 
scattering event on the iso-delay surface, while the remain- 
ing fraction was scattered by a number of interactions with 
the dust grains along the hght path I within the cloud sys- 
tem. The path length I is defined as the cumulated distance 
traveled from the central source through the single scatter- 
ing point and to the outer boundary of the dust cloud in 
the direction of the observer across regions with n 7^ (see 
Fig-EJ. In the most general case, / depends on the ip, z, t co- 
ordinates of the scattering point (see Fig.Qfor the meaning 
of ip). The inclusion of attenuation is achieved introducing 
a exp [— Ta(<^, 2, t)] factor in the definition of f{t): 



f(t) = Cext C 



n{(p,z,t)e dip dz (11) 



-ct/2 Jo 

where the attenuation optical depth Ta is defined as 



Ta{(p,Z,t) = C, 



ext I n{l) \dl\ 



(12) 



For one of the mean value theorems (see for example 
iRade fc Westergreen|[l999l) . one can also write: 



Cext LO C _ 

fW = 7; e 



G{t) 



where: 



J-ct/2J0 r ' ' 

I-Z2 ^ '*('^' ^' '^'^ 



dip dz 



(13) 



(14) 



and G{t) is a time dependent function which includes 
the properties of dust and dust distribution. From Ea. ll3l it 
is clear that, with the introduction of attenuation, the over- 
all LE luminosity increases for growing density values until 
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the optical depth reaches a critical limit. After that, self- 
absorption prevails and any further density enhancement 
causes the luminosity to decrease. 

In general, r^ff is time and geometry dependent and 
it can be regarded as a weighted optical depth for the LE 
at any given time. For example, it is easy to demonstrate 
that for a SN immersed in a perpendicular dust slab with 
uniform density, Ta does not depend on tp and z coordi- 
nates, since I = D + ct. Therefore the solution of Eq. 1141 
becomes straightforward and it gives Teff = Td + Cext no ct, 
which implies that the deviation from the single scattering 
approximation grows with time. A few more realistic exam- 
ples of Teff calculations for the V band are presented in 
Fig. |51 where I have considered four different geometries: a 
uniform spherical shell extended between R and R + AR, a 
distant perpendicular sheet placed at distance D and with a 
thickness A_R, a distant spherical cloud with a radius R and 
located at a distance D in front of the SN and a spherically 
symmetric wind with n(r) = no (Ri/r)^, which extends be- 
tween Ri and R2 > Ri- The first case can be considered as 
an approx imation to the so-called contact discontinuity (see 
ISugermanlF 2003V This can be produced by the mass loss un- 
dergone by the la progenitor binary companion, through the 
interaction between the slow dense red supergiant wind an d 
the fast blue supergiant wind jEmmering fc Chevalieiil989ll . 
For the thin shell and perpendicular sheet I have adopted 
values which are very similar to those of Sugerman (2003; 
see his Table 4, second B case and first / case for SNe re- 
spectively), for the spherical cloud I have use d a typical siz e 
of molecular clouds systems (see for example iMathi ^[200011 . 
while for the stellar wind I have allowed for an inner dust- 
free cav ity with rad ius 7?i=0.5 lyrs (~5xl0^^ cm, see for 
example 'Dwok"l983^ and an outer radius R2/Ri = 15. As for 
the dust parameter s, I have adopt ed the values r eported by 
IWeingartner fc Drai nc ( 2001) and iDraind for the V 

passband (see also Table I^J. 

For all geometries, the number densities were set in or- 
der to produce the same optical depths Td, i.e. the same SN 
extinction. The relevant parameters are summarized in Ta- 
ble As one can notice, for the stellar wind case the higher 
optical depths imply dust densities which are rather unreal- 
istic for this geometry and we re included for completeness, 
also following IChevaheil lll986ll . 

As far as the relation between Tef / and is concerned, 
it must be noticed that this depends very much on the dust 
geometry. For example, for a system with no dust along the 
line of sight, Td—O, while Teff can assume any value. This 
clearly shows that, in general, the optical depth of a LE 
cannot be expressed in terms of the optical depth of the SN 
as seen from the observer. This is approximately true when 
the density is low or for some particular geometries only. For 
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the case of an isolated and distant cloud), the generation of 
random photons on the whole unit sphere would be not only 
a waste of computing time, but would also be inconsistent 
with the forced scattering technique, which I will introduce 
in the next sub-section. This problem can be solved emit- 
ting photons only within the solid angle Qe subtended by 
the cloud at the SN position. The fact that photons are ar- 
tificially generated within a cone and not on the whole unit 
sphere is taken into account during the normalization of the 
light curve (see Sec. 16.41 Ea. llSH . 
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Figure 5. Example r^f j calculations for T|J=0.03 (solid line), 1.0 
(dashed line) and 2.5 (dashed-dotted line) in SSA approximation 
and V passband. Geometrical parameters are indicated in each 
panel. For comparison, filled symbols trace the corresponding MC 
solutions (see Sec.|7J. 



example, for the slab case seen before, r^ff — Td for ct <^ D. 
As we will see, Te// actually depends not only on time and 
geometry, but also on the dust properties. 



6 MONTE CARLO SIMULATIONS FOR 
MULTIPLE SCATTERING 

The idea behind Monte Carlo simulations is to follow each 
single photon packet along its random walk until it finally 
escapes the dust cloud. This procedure implicitly includes 
the treatment of multiple scattering and it is therefore suit- 
able for all ca ses where the du st optical depth is large. 

Following 'chevalier' (l98^, here I will adopt the proce- 
dure outlined bv .Witt l,1977ll . recalling the fundamentals of 
the method and describing in detail the aspects which are 
more specific to our case. 

6.1 Isotropic and beamed emission 

The whole process starts with the photon emission. Due to 
the dimensions of the dust clouds I am going to consider, I 
can assume that the SN is a point source lying at the origin 
of a coordinate system x,y,z. If 5Ri is a random number 
uniformly generated in the range 5Ri 5^ 1, the direction 
cosines are given by: 



Wo 
Uo 

vo 



2K, - 1 



cos [7r(25R,:-| 



v/l-^o sin[7r(2K>+i - 1)] 

and o bey to the usual condition uq + Vq -j- Wo = 1 (see 
IWittligf^ . Eqs. 5, 6 and 7). When the scattering material, as 
seen from the central source, shows significant voids (like in 



6.2 Forced first scattering 

To avoid wasting time with photons directly escaping the 
cloud, I have followed the common pr actice using the so 
called forced first scattering The basic idea be- 

hind this technique is to scatter each emitted photon, while 
the direct escape of non interacting photons is then taken 
into account using a weighting factor Wo = 1 — exp~^i . 

Once the random direction is generated, one needs to 
compute the total optical depth along the light path r: 



ri 



n{r) Cext dr 



where rmax is the maximum dimension of the cloud or 
system of clouds along the emission direction. The next step 
is the random generation of r ^ n at which the first scatter- 
ing will occur. This is achieved using as the probability 
density function, applying the inverse method and solving 
for T. This gives the following expression: 



r = --ln[l-3ff(l- 



which corresponds to Eq. 14 of lWitli 11197711 . except from 
the fact that the latter has the incorrect sign. 

Once this is done, one needs to determine the spatial 
position r where the scattering occurs. For simple dust ge- 
ometries this can be obtained inverting the analytical re- 
lation which links r to the particle density function, while 
more general cases can be solved num erically (see for ex- 
ample |Fisc^^^ennin^^^%r3 ^9^, at the expense of 
computing time. 

Irrespective of the adopted method, one has then to 
compute the photon direction after the scattering event, in 
order to be able to decide about its immediate destiny. For 
this purpose, as in the case of the analytical solutions, I have 
adopted the H-G scattering phase function. Again, using the 
inverse method, and solving for /ii one finally obtains: 



Ml = 



(1+ff 



(l + 2ffK-<?) = 



(15) 



19771) , except from 
'^^ term was lost. 



which corresponds to Eq. 19 of lWit 
the fact that the exponent 2 on the {1 — g 
For g — 0, Eg. llSl has to be replaced by ni = 25R — 1. 

Once /ii is generated according to Eg. 1151 the new direc- 
tion cosines {ui,vi,wi) need to be computed. This is done 
by the following coordinate transform ation (see for example 
iBianchi. Ferrara fc Giovanardilll996l Appendix A4): 



"0 



Ml = Uo cos t> 

vi — Vo COS 9 
till = Wo cos 9 — -^1 — Wq 



+ 
+ 



sing cos ( 

VQ-WQ 



■■ sm o cos ( 



sm o cos ( 
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where 9 is the scattering angle defined by ni = cos 9, 
and (j) is the scattering azimuth, from which I assume the 
scattering efficiency to be independent^. This is generated 
simply as (j) = n (25R— 1), which gives a uniform distribution 
in the range — vr ^ ^ tt. 



6.3 Successive scatterings and final photon escape 

The optical path to the next possible scattering can be com- 
puted with considerations which are similar to those I made 
for the first scattering, except from the fact that now pho- 
tons are not forced to interact with the dust. Under this 
condition, the optical depth at which the next scattering 
event will occur is simply given by r = — In 5R. Whether the 
photon will escape or undergo further scatterings can be de- 
cided comparing r with T2, which is now computed along 
the scattered direction from the scattering point to rmax- If 
r > T2 the photon will escape and, to compensate for the 
first forced scattering, a weight W = uj Wo will be given to 
it. The dust albedo uo <1 accounts for the fact that during 
the scattering some photons can actually be absorbed. 

In the case oir the photon packet will be scattered 
again and hence the whole process has to be repeated, start- 
ing with the calculation of the scattering point, scattering 
direction and so on, until it finally escapes. In the general 
case of A*' ^ 1 scattering events, the final weight will be 



that the total optical path \s L — I + "Ylii^i Having this 
in mind, the arrival time ta can be easily computed as: 



W = tj™ Wo 



(16) 



Once the emerging fraction of a given photon packet 
escapes, there are a series of parameters that are relevant 
for its temporal and spatial classification: a) the position 
Pn{xn ,yN , zn) of last scattering; b) the direction cosines 
after last scattering (ujv, vn, wn)\ c) the time delay accumu- 
lated by the scattered photons with respect to those which 
were emitted at the same time but traveled directly to the 
observer. All these numbers are needed to reconstruct inte- 
grated luminosity and surface brightness image (or profile) 
at a given time and from a given viewing direction. 

While parameters a) and b) are immediately available, 
the calculation of time delay requires a little bit of more 
discussion. 



6.4 Light curve reconstruction 

Escaping photons can be seen only by an observer who lies 
on the line of sight coinciding with the last scattering direc- 
tion p{uN ,vn ,wn)- If d is the vector which goes from the 
coordinate system origin to the observer, r{xN,yN, zn) de- 
fines the position of last scattering Pjv and I = Ip is the vec- 
tor which goes from Pjv to the observer, one can clearly write 
d = r + l. Then if I set k — r op = {xnUn + yNVN + znWn), 
I obtain = + 2kl + , which can be solved for /. One of 
the two solutions of this equation is not physical and there- 
fore the light path from the last scattering to the observer 
is I = —k + y/k'^ + — r^, where = x% + y% + z%. For 
d^ r this gives I ~ d — k. Taking into account the random 
walk which was traveled before the last scattering, one finds 



^ Note that the equivalent Eqs. 22 given bv IWitj ^1977^ are af- 
fected by several typos. 



JV 

Eri I — d 

— + 
c c 



JV 

1=1 



(17) 



where is the emission time. This has to be generated 
according to the input SN light curve L{t), which in this case 
can be regarded as a probability density function. Applying 
the inverse method I can write: 



5R = 



1 



Lq Atsf 



L{t) dt 



Since, in principle, one does not have an analytical de- 
scription of L(t), the integral has to be evaluated numeri- 
cally for a series of te values and the result can be easily 
interpolated and inverted, so that entering the uniform ran- 
dom variable 5R one gets te with the correct distribution. Of 
course, in the case of flash approximation (see Eq. 0, the 
emission time is simply generated as tg = AtsN 

At this point, the light curve can be progressively com- 
puted just accumulating the weighted photon packets in 
time delay bins AT along a given direction and within a 
given solid angle Afl. Now, if p{uN ,vn ,wn) is the pho- 
ton escaping direction and s{sx, Sy, Sz) is the versor which 
identifies the line of sight, the angle tp formed by the two 
vectors is defined by costp = p o s. Photons are counted 
if cosi/i ^cosi/);, where the limiting angle can be defined as 
cos-i/); = 1 — q and AQ — 2nq, with ^ g ^ 1. 

The bin size AT, which in the end will give the time 
resolution, depends on the total number of generated packets 
and the signal-to-noise ratio one wants to achieve. The same 
applies to AQ, which will govern the angular resolution. 

Finally, in order to follow the same procedure I adopted 
for the analytical solutions, the resulting light curve needs 
to be normalized to the SN maximum light luminosity. This 
is achieved dividing it by the number of photons Fo the 
observer would receive by the SN at maximum, within the 
solid angle ASl and in the time interval AT, which is clearly 
given by 

To = iV,,„t ^ (18) 

AlSJV lie 

For spherically symmetric models, escaping photons can 
be both generated and collected all over the unit sphere 
(An = fie = 47r), and hence the normalization is obtained 
dividing the photon counts by Nphot AT/ AtsN ■ 

6.5 Echo surface brightness reconstruction 

Escaping direction cosines and last scattering position can 
be used to derive the spatial energy distribution projected 
onto the plane perpendicular to the line of sight, which co- 
incides with the plane of the sky. This can be achieved sim- 
ply transforming the coordinates of the last scattering point 
Pn(xn ,yN , zjv) in the new system 0[X, Y, Z), which is ob- 
tained rotating by an angle (5 around the z-axis and by an 
angle 7 around the transformed y-axis. The rotation angles 
P and 7 are defined by the escaping direction cosines as: 



1/3 = vo/\Jl — vol; cosP = uo/\/l 



(19) 



sm7 = Wo; 



57 = yjl-wl 
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In this reference system, the plane of the sky coincides 
with the Y, Z plane, while the X-axis coincides with the line 
of sight. The transformation equations are as follows: 

X = cos a cos /3 + j/jv sin a cos /3 + zjv sin (3 
Y = —xn sin /3 + j/jv cos /3 

Z = —xn cos I3 sin 7 — j/jv sin /3 sin 7 + zjv cos 7 

Of course, contrarily to what happens for the numeri- 
cal integration of Eq. |S| MC simulations do not allow direct 
calculation of the echo surface brightness at a given time 
t, since the arrival times ta are scattered across the whole 
range allowed by the dust geometry. For this reason, one 
needs to perform a selection on the outcoming photons and 
collect them only if their arrival times satisfy the condition 
t - AT/2 < ta < t + AT/2. Due to the large difference be- 
tween the overall echo duration and the time resolution AT, 
this means that only a small fraction of the total number of 
generated photons will fall in the required time range. 

In order to achieve a reasonable signal to noise ratio 
in the final results, this implies that a very large number 
(10^ — 10^) of photons must be generated. 



7 RESULTS OF MC SIMULATIONS 

I have implemented the concepts discussed in the previous 
sections in a MC code, which I have used to study the ef- 
fects of multiple scattering on the outcoming LEs for the 
four dust geometries described in Sec. |3 In all calculations 
I have adopted Cext=5.21xl0~'^^ cm^, p=0.54 and uj^OM, 
which are typical of a cano nical Ry =3.1 Milky Way mix- 
ture in the V passband (see lDrainell200 j and Sec. |H| here). 
For non spherically symmetric geometries, I have used Ar2= 
2xl0~^7r {q=10~^), which corresponds to a collecting beam 
with an angular semi-amplitude of ~2.6 degrees. 

The basic effects of multiple scattering are illustrated 
in Fig. I^for different values of Td- While the MC simulated 
luminosity is only a few percents fainter than the single scat- 
tering solution in the low optical depth case (rd=0.03), this 
difference becomes more and more marked as the dust den- 
sity grows. For Td=0.1 the deviation is ~0.1 mag and it 
reaches about 1 mag for Td=1. 0. This is in very good agree- 
ment with the results found bv lCheva licr ( 1986) in the case 
of type II SNe surrounded by dust distributed with a 
density law. 

This is clearly due to the growth in the number of mul- 
tiple scatterings. In fact, while the number of single scatter- 
ings tends to increase the echo luminosity, further scatter- 
ings contribute to the auto-absorption and therefore work 
against any brightness increase. As a consequence of these 
two competing mechanisms, the echo luminosity keeps grow- 
ing for increasing values of n until the fraction of photons 
which undergo pure single scatterings reaches 0.5. The scat- 
tered flux increases more and more slowly as the dust optical 
depth approaches some critical value (t^ —1), after which 
the fraction of photons with Nscat ^ 2 becomes larger than 
0.5 and the overall luminosity decreases for any further in- 
crease of Td- The effect of a further dust density rise is illus- 
trated by the gray circles, which were obtained for Td=2.5. 
The exact details depend on the specific dust distribution 
but, for example, in the case of the thin shell, under these 
circumstances Nscat~l in ~15% of the cases only and this 



turns into an overall echo luminosity which is more than 1.5 
mag fainter than the corresponding single scattering solu- 
tion. Very similar results are obtained for the distant sheet 
and the cases, while for the spherical cloud the be- 
haviour is different, due to the fact that the intersection 
between the LE paraboloid and the cloud tends to rapidly 
decrease. 

While this is not foreseen in the plain single scattering 
approximation, this is qualitatively predicted by the SSA 
approximation (solid curves), which gives reasonably good 
results for r^ff < 1. In general, the simulations show that 
SSA successfully reproduces the early phases under a wide 
range of geometries, while it progressively deviates from the 
MC solutions as time goes by, and the deviation is larger 
for higher values of r^ff. The light curves produced by SSA 
are systematically fainter than multiple scattering solutions. 
This is a consequence of the assumption on which SSA is 
based, i.e. that attenuation acts as a pure absorption, which 
is approximately true at low optical depths only. In fact, 
when the density grows, photons which are deviated by the 
single scattering trajectory still have a probability of reach- 
ing the observer through further scatterings and can, there- 
fore, contr ibute to the LE a t later times. This effect, already 
noticed bv lChevalierl (Il98dl , translates into a light curve flat- 
tening and it can be seen as the consequence of an increased 
average time delay accumulated by photons. The net effect 
is to reduce the effective LE optical depth with respect to 
Tsff implied by the SSA approximation (see Fig. |K1 filled 
symbols). The simulations show that this is more effective 
when the dust is confined closer to the SN and, for a given 
Teff, when the dust density is higher. 

More subtle is the outcome of multiple scattering on 
the radial surface brightness profile. For a uniform dust dis- 
tribution, the simulations show that the resolved LE would 
always appear as a disk with a small surface brightness in- 
crease going outward and the gradient tends to become more 
pronounced for larger values of the dust optical depth. How- 
ever, the effect is practically negligible. 

Finally, the dust optical depth causes a progressive 
breaking of the relation between r and p expressed by Eq. |5] 
which, in principle, allows the direct calculation of the dust 
distance from the SN from the resolved LE image. In fact, in 
the high density regime the contribution by multiple scatter- 
ings tends to fill all the region internal to the iso-delay sur- 
face, mimicking the presence of dust close to the SN itself. I 
must notice, however, that the effect is not really pronounced 
and in most cases it would not be observable, the multiple 
scattering contribution being lost in the background noise. 



8 LIGHT ECHO SPECTRUM 

So far I have restricted the discussion to one single wave- 
length. But, in principle, a LE can be studied at any wave- 
length, provided that the dependencies of the relevant pa- 
rameters, namely Cext, lo and g, are known. If this is the 
case, the light echo spectrum can be computed using Eq. Q 
where L{t) is replaced by the spectrum S{X,t). 

For this purpose I ha v e ado pted the resu l ts pub lished 
by IWeingartner fc Drainel J200J) and iDraind J2003l) for a 
Milky Way mixture of carbonaceous and silicate grains with 
a canonical Rv=S.l. The behaviour of the relevant param- 
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Figure 6. MC solutions for four different dust geometries for r£j=0.03 (triangles), 0.1 (squares), 1.0 (circles) and 2.5 (gray circles). For 
comparison, single scattering and SSA solutions are also indicated (dotted and solid lines respectively). For clarity, the SSA solution for 
T;j=2.5 has been traced with a thick gray line. 



eters in the optical spectral range is shown in Fig. where, 
for comparison, I have also plotted the functions computed 
for higher values of Ry- As one can see, the dust albedo 
changes very little across the optical range (0.60^ u) ^0.67), 
while the forward scattering degree shows a more significant 
variation. However, the strongest wavelength dependency is 
shown by the extinction cross-section whose behaviour, in 
the wavelength range 4000-9000 A, is well approximated by 
a A"" law, with a=1.35. 



These functions were implemented in the codes by 
means of polynomial interpolations to the tabulated values''. 

8.1 Single scattering 

In single scattering approximation, as in the case of the hght 
curve, the echo spectrum shape is independent of no, which 

^ Tabulated functions can be found at the following URL: 
: / /www . astro ■ princeton ■ edTi/~dralne/dust /dustmix ■ html | 



|http 
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Figure 7. Behaviour of albedo (lower panel) forward scattering 
degree (middle panel) and extinction cross section (upper panel) 
for a Milky Way size dust distribution (Weingartner & Drainc 
|2001, Drainc 2003) and different values of Ry. Filled dots indicate 
the values at the effective wavelengths of B and V passbands. 



governs the global luminosity only, while dust geometry and 
dust properties effects are included through the impulse re- 
sponse function f{t). Now, the SN luminosity decays very 
fast and thus the bulk of optical radiation is emitted in a 
very short time range'*, which I will indicate as Ati,. As a 
consequence, what really counts for the shape of the output 
spectrum observed at time t is the behaviour of / between 
t — Ati, and t. Therefore, the geometrical dependency is ex- 
pected to be generally very mild ai t ^ Atb- Exceptions to 
this might be encountered, for example, during the onset of 
a real echo, in which / shows a pronounced transient phase. 
For the same reason, also variations on the forward scatter- 
ing degree g are expected to produce significant changes in 
the overall flux, but not in the spectral shape. 

This is confirmed by the numerical solutions, for which 
I have used the s pectroscopic data set of SN 1998aq 
jBranch et al J 120031) . This includes 29 spectra obtained in 
the phase range —9'* 5; i ^ +241'* and covers the wave- 
length range 3720-7160 A. An example of LE spectrum is 
shown in Fig. |S] it is immediately clear that it does not re- 
semble the SN spectrum at maximum, even if it has roughly 
the same colour (see also next section). For comparison, in 
Fig. |H] I have also plotted the time integrated spectrum of 
SN 1998aq (gray line). This is clearly redder than the actual 
LE spectrum, mainly due to the wavelength dependency of 

Cext- 

As I have already pointed out, in single scattering ap- 
proximation an increase of the dust density produces a pro- 
portional increase in the overall echo flux, but no variation in 

* For a normal la, 90% of the radiation between 3700 A and 
7200 A is released in about 0.3 yrs. 



Figure 8. Light echo spectrum computed at i=2 yrs in single 
scattering approximation for a thin shell with i?=20 lyr, AR=2 
lyr (black line) and fa Tgjj=0.03 (rao=25 cm~^). For compar- 
ison, the solid gray line and the dotted line trace the time inte- 
grated spectrum and the maximum light spectrum of SN 1998aq 
(Branch et al. 2003), respectively. For presentation all spectra 
have been normalized to their integrated luminosity in the wave- 
length range 3720-7160 A. 

its spectral appearance. In order to investigate optical depth 
effects, one in fact needs to use the multiple scattering ap- 
proach. This is of course possible, even though the number of 
photons required to reach a reasonable signal to noise ratio 
in the output spectrum becomes very large. Practically this 
allows echo spectral synthesis in spherical symmetry only 
and, therefore, geometrical effects must be evaluated using 
broad band colours instead (see next section). 

8.2 Multiple scattering 

The procedure to compute a multiple scattering spectrum 
is identical to the one I have described for the light curve, 
with the exception that now the emission time te and the 
emitted wavelength Ae have to be generated according to 
the input spectral distribution S{X,t). This can be easily 
achieved using the inverse theorem. Once the photon has 
completed its random walk through the dust, it is finally 
classified according to the arrival time ta and counted in 
the proper wavelength bin, whose amplitude AA fixes the 
spectral resolution. 

The results for a thin sheU (i?=20 lyr, AR=2 lyr) and 
t—2 yrs are shown in Fig. |5| for three different values of 
Td- As expected, the introduction of multiple scattering pro- 
duces significant reddening in the spectrum. As a matter 
of fact, the spectra calculated for higher densities can be 
obtained directly from the one computed for the lowest Td, 
applying the same wavelength dependency of Cext used for 
the LE calculation and a suitable optical depth re//. An 
example is presented in Fig. |5| (solid thin line), which was 
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Figure 9. Synthetic LE spectra for a thin spherical shell (i?=20 
lyr, AR=2 lyr) at t=2 yrs for increasing values of (dark to 
light thick lines). The thin line is the low optical depth spectrum 
reddened with a suitable extinction (see text); for clarity the spec- 
trum was slightly shifted downward. For comparison, the dotted 
line traces a spectrum computed for t^=2.5 with the SSA ap- 
proximation. For presentation all spectra have been normalized to 
their integrated luminosity in the wavelength range 3720-7160 A. 



obtained using tJ^j. = 1.95, i.e. the value derived from multi- 
ple scattering light curve model (see Fig.l^. For comparison, 
in Fig. have also plotted the spectrum obtained in SSA 
approximation for rj = 2.5 (dotted line). This appears to 
be too red with respect to the corresponding MC solution 
and the reason is that SSA tends to overestimate the LE 
optical depth by amounts which become larger at shorter 
wavelengths (see also next section). 



9 LIGHT ECHO COLOUR 

An alternative to spectral synthesis is the much less expen- 
sive calculation of broad band light curves, which can be 
combined to produce synthetic colours. In order to explore 
the potentialities of this method, I have studied the case of 
the {B — V) colour, since for a la the LE is clearly brighter 
in these passbands than at redder wavelengths (see Fig. (HJ. 



9.1 Single scattering plus attenuation 

A first analytical estimate can be derived in SSA approx- 
imation, assuming that the input light curve is a flash. 
According to what I have shown in Sees. |21 and |3 un- 
der this assumption the echo flux at a given wavelength is 
F = Lq no Cext ^tsN e~^^ff G{t), where G{t) is a time and 
wavelength dependent function related to geometry and dust 
properties. Applying this to the case of B and V passbands, 
and indicating with (B — V)^ax the unreddened colour of 
the SN at maximum, one can estimate the LE colour as: 



Table 2. Basic parameters for B, V light curves and dust mix- 
tures. Value s for Ce x t, ^ a nd g are from Wcingartncr & Drain^ 
i200ll) and iDraind <2003l) for Rv=3.1 (first two rows) and 
Ry=5.5 (last two rows). Extinction cross section Cext is given 
per Hydrogen atom. 
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0.76 


0.60 





(B-V) 



(B - VfLax ~ 2.5 log 



l-SN 



GYxt ' 



;^ At;; 



+ 



1.086 {T^ff 



'eff 



(20) 



where, due to the mild wavelength dependency of g (see 
Table I have assumed G^(t) ^ G^{t). For low values of 
Teff, after substituting the relevant parameters in Eg. 1201 
one gets {B - V) ^ {B - V f^ax ^ -0.1. Therefore, when 
the single scattering approximation is valid, the LE colour 
is roughly time and geometry independent. More generally, 
when the LE optical depth is not negligible, Eo. 1201 implies 
that the resulting colour is reddened by auto-absorption 
{r^ff > T^ff)- Moreover, it is directly related to the LE 
optical depth. In fact, using the values reported in Table |51 
Eg. l2Ul can be rewritten as: 



{B-V)^{B 



V)° 

' J max 



+ 0.36 r^ff 



(21) 



Thus, the observed {B — V) colour can be used to esti- 
mate the effective optical depth of a LE, so that this is not 
a free parameter. Actually, numerical simulations show that 
the ratio G^{t)/G^{t) can vary with time and geometries 
between 0.9 and 1.1, so that the previous simplifled equa- 
tion is expected to produce maximum errors of about ±0.1 
on the derived {B — V) colour. 



9.2 Multiple scattering 

As we have seen in Sec. |7| when multiple scattering is at 
work, the echo luminosity does not scale linearly with the 
density and therefore the observed colour is expected to 
show a dependency on the dust optical depth. In partic- 
ular, for increasing values of particle density, light curves 
appear to get fainter and fainter with respect to the corre- 
sponding single scattering solutions (see for example Fig.jSJ. 
Since this effect is supposed to get stronger for higher values 
of Cext (i.e. at shorter wavelengths), one expects the echo 
colour to become redder and redder as the density grows, as 
qualitatively foreseen also in the SSA approximation. 

To explore these ideas in a more quantitative way, I 
have run a series of calculations using the parameters listed 
m TableU which were derived from the template light curves 
presented in Fig. |21 and from Draine's models. The results 
for the thin shell and the r^^ wind are shown in Fig. 1101 for 
three different values of , while Fig. llll presents the cases 
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Figure 10. Multiple scattering (B — V) colour curves for the thin 
shell (circles) and wind (squares) for three different values 
of T|j. The corresponding SSA solutions are traced by solid and 
dashed curves respectively. 



Figure 11. Multiple scattering (B — V) colour curves for the per- 
pendicular sheet (circles) and spherical cloud (squares) for three 
different values of t^. The corresponding SSA solutions are traced 
by solid and dashed curves respectively. 



of the distant sheet and the spherical cloud. For low optical 
depths the MC solutions are practically identical to the SSA 
solutions (solid and dashed curves) while, for larger densi- 
ties, the calculated colour curves turn globally redder and 
redder for both geometries and progressively deviate from 
the SSA solutions. In general, these latter are systemati- 
cally redder than the multiple scattering solutions and they 
give reasonable results only for r^ff <1. 

As far as the time dependence is concerned, the degree 
of forward scattering plays an important role. In fact, the 
random walk of a photon through the medium is influenced 
by the scattering phase function, in the sense that for high g 
values, a photon tends to follow its emission direction, since 
small scattering angles are favored, and therefore its random 
walk will be shorter than for lower g values. This means that 
the light curve is expected to declines faster at bluer wave- 
lengths, where g is higher (see Fig. |7|l. On the other hand, 
since Cext also increases bluewards, this produces the oppo- 
site effect, implying a light curve flattening, which is more 
efficient in the blue. As the si mulations show, for the Milky 
Way mixture we have adopted iWeingartner fc Drainel200ll 
iDraindboOSl) the two mechanisms tend almost to compen- 
sate each other, giving a rather flat time evolution, mostly 
driven by geometrical effects. As a matter of fact, the exact 
behaviour turns out to be pretty sensible to the wavelength 
dependency of g. Both SSA and MC calculation show that if 
one assumes dg/dX=0, the {B — V) colour tends to become 
bluer with time. Of course, the time evolution will also de- 
pend on the dust distribution, mainly through the allowed 
incoming angle range variation and the changes in Te//. For 
example, when the dust is far from the SN, as in the case of 
the perpendicular sheet (Fig. ^3 circles) , this does not vary 
significantly with time, and therefore the colour curve is ex- 
pected to be flatter than for dust close to the central source. 



This is also true for the distant spherical cloud (Fig. 1111 
squares), but the fact that Te// decreases as time goes by 
causes the LE to become progressively bluer. 

The main conclusion is that the LE colour is practically 
time and geometry independent, especially after the initial 
onset phase {t >2 yrs), provided that the dust density is 
uniform. Therefore, any significant change in the observed 
colour has to be interpreted as a variation of Te// which, in 
turn, would signal the presence of density fluctuations. 

As in the case of SSA approximation, one can use the 
approximate expression given by Eg. l21l to estimate the ef- 
fective multiple scattering optical depth directly from the 
observed {B — V) colour. It must be noticed that the as- 
sumption that T^ff/T^ff = Cfxt/CYxt^ which holds in SSA 
approximation, is not strictly true when multiple scattering 
is at work. In general, in fact, this depends on the particular 
geometry and tends to be slightly larger than the cross sec- 
tions ratio. Nevertheless, the simulations show that Eq. 1211 
gives reasonable results for r^jf ^ 2, where now e'^'ff is 
defined as the ratio between the multiple scattering solution 
and the corresponding single scattering solution. 



10 EFFECTS OF A DUST MIXTURE CHANGE 

In all calculations, I have so far adopted a canonical Rv^S.l 
Milky Way dust composition. To illustrate the effects of a 
change in the dust mixture, I have run a series of calculations 
using Draine's model corresponding to Rv~5.5. As shown in 
Fig.[7| the three dust parameters increase at all wavelengths. 
In the single scattering approximation, the consequences on 
B and V LE luminosities and colours can be anticipated 
from the following considerations. Neglecting the small vari- 
ation in the forward scattering degree (see below), the LE 
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Figure 12. Multiple scattering {B — V) colour curves for the 
thin shell case and for Ry=5.5 (symbols). The densities have 
been reduced in order to give the same optical depths as in 
Fig. ESI For comparison, the corresponding solutions for Ry=3.1 
are also plotted (lines). 



luminosity is expected to change proportionally to uj C'ext 
(see Eq.|2l. Using the suitable values (see Tablel^J this turns 
into an increase of 0.29 and 0.44 mags in B and V respec- 
tively, while the {B — V) colour becomes 0.15 mag redder. 

When multiple scattering becomes important, the exact 
outcome is not so obvious. As I have shown in Sec. Q a 
dust optical depth enhancement produces an increase in the 
LE luminosity until this reaches some maximum value, after 
which absorption prevails on scattering and the luminosity 
starts to decrease again. Therefore, increasing the scattering 
cross section causes an earlier luminosity saturation for any 
given dust density. 

On the other hand, an albedo increase translates into 
smaller losses due to absorption, and therefore into a 
brighter light reverberation. The same mechanism tends also 
to delay the luminosity saturation, since higher values of the 
albedo turn into a smaller efficiency of multiple scattering in 
absorbing the radiation. The effect is going to be more pro- 
nounced at higher densities since, in those conditions, a fair 
fraction of the photon packets which reach the observer have 
undergone multiple scattering and therefore their weight, 
given by Eg. 1161 strongly depends on the albedo. Finally, an 
higher forward scattering degree tends to increase the lumi- 
nosity, especially at early phases, and to enhance the time 
dependency (see Sec.|HJ. 

To illustrate the combined action of these mechanisms, 
in Fig. ll2I I have plotted the comparison between the (B — V) 
colours obtained for Ry = 3.1 (lines) and Ry = 5.5 (sym- 
bols), for the thin shell case (_R=20 lyrs, AR=2 lyrs). For the 
high Rv solutions, the densities no were adjusted in order to 
produce the same optical depths tJ , i.e. 0.03, 1.0 and 2.5. 
As expected, for the low density case, the Rv=5.5 solution 
(triangles) is about 0.15 mags redder than the i?v'=3.3 one 



(solid line), the difference being larger at earlier phases due 
to the forward scattering effect. Then, at higher densities, 
the two solutions (squares and dashed line) become similar, 
implying that the optical depth effect is weakened by the 
albedo effect. Finally, for rj'=2.5 (circles and dotted line) 
the -Rv=5.5 solution is bluer than the Rv=3.1 one. As a 
consequence, it is not possible to deduce the optical depth 
and the extinction law from the observed LE spectrum, even 
in the hypothesis the unreddened SN spectra are available. 
At least one of the two has to be assumed in order to derive 
the other. In fact, at least for the two different dust mixtures 
we have used here {Rv=3.1 and 5.5), one can produce very 
similar results just using different values for the LE optical 
depth. 



11 LIGHT ECHO POLARIZATION 

Due to dust scattering, resolved LEs are expected to show 
strong linear polarization, which should reach its maximum 
degree for the dust which produces a scattering angle 9=tt/2 
JSparkJl994 1997). In single scattering approximation, this 
occurs for dust lying on the SN plane, at a distance ct from 
the SN itself. 

In this section, I want to address the effect of multi- 
ple scattering on the observed polarization. A similar prob- 
lem, even though in rather different contexts, has been 
approached using MC te c hniques by several authors (see 
for example lHoflichlll99ll. iFischer. Henniiig fc YorkdITggl. 
ICode fc WhitnevI Il995l. IWood et alJ Il996l Bianchi et al. 
1996), and here I basically follow the same procedure they 
have adopted. 

The fundamental concepts are those exposed by Chan- 
drasekhar (1950, Chap. I; hereafter C50). In this view, a 
photon packet is characterized by a Stokes vector S = 
{I ,Q,U,V), which at each scattering event undergoes a 
transformation described by the following equation (C50, 
Eq. 210): 

S ^ L(tv -ii) ■ R - L{-ii) ■ So (22) 

where So denotes the incoming Stokes vector, R is the 
so-called scat tering matrix and L is a Mueller rotation ma- 
trix (see e.g. iTinbergeiil ll996L Chap. 4) which transforms 
the Stokes vector from one reference system to the other. 
Following C50, I use as reference plane that defined by 
the photon direction and the 2-axis, so that the required 
transformations are a counter-clock rotation of amplitude 
ii , to bring the Stokes components into the scattering plane, 
and a clockwise rotation of amplitude vr — 12 to recompute 
the Stokes co mponents with the new reference plane (see 
I Van De Huis^ lTgSOy 

The explicit terms f or L can be derived fr om C50, 
Eq. 186 (see for example ICode fc WhitnevI Il995l Eqs. 3), 
while approximate expressions for R term s in th e case of 
H-G phase function can be found in |WhitJ(|l97i). 

11.1 Single scattering polarization 

The effect of a single scattering on an unpolarized photon 
packet (Qo = (7o = Vb = 0) emitted by the central source, 
can be immediate ly computed using the approximate ex- 
pressions found bv lwhitg ^1979^ and this gives 
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Pi ^pi 



1 + cos2 e 



(23) 



0.6 



where pi, the peak Unear polarization, is of the or- 
der of 0.5. This result, first applied to LEs by Sparks 
(1994, Eq. 10), implies that the polarization is maximum 
for 8=tt/2, i.e. on the SN plane. 

As for the polarization angle, it is easy to verify that the 
scattering pla ne and the elec tric field are perpendicular. As 
pointed out bv lSpark^ jigged) , this implies that the polariza- 
tion direction is, at any given position, tangent to the circle 
centered on the SN, and this can be used as a distinctive 
feature to identify a LE. 

Using Eq. |7| in Eq. 1231 one can easily compute the ex- 
pected radial polarization profile (see Fig. 1131 solid line), 
which can be used to compare with the multiple scatter- 
ing solutions I am going to present in the next sub-section. I 
note that in single scattering approximation the polarization 
degree is independent from no, g, Cext and uo. 



11.2 Depolarizing effect of multiple scattering 

The implementation of polarization calculations in a MC 
code is fairly simple and it does not differ substantially from 
the one I have used to compute the echo surface brightness. 
After generating an unpolarized photon packet, its polariza- 
tion state is followed along its random walk through Eas. l22l 
In order to take into account the forced nature of first scat- 
tering, the initial Stokes vector is set to So = (Wo, 0,0,0), 

where Wo = 1 — exp{—Ti) (see Se c. 16.211. 

As usual in MC methods (see lKattawar fc PlasjlT96d) . 
the post-scattering Stokes parameters are weighted by the 
distribution which they were sampled from, which in this 
case is Pi itself, i.e. the H-G phase function. In our case, 
all Pi (i=l,2,3,4) terms of the scatterin g matrix contain Pi 
as multiplicative factor (see lWhiteligTgl) and, therefore, this 
corresponds to setting Pi = l. Finally, in order to account for 
the albedo effect and to make the radiation transfer consis- 
tent with the one I have used for the luminosity calculations, 
one would need to normalize the new Stokes parameters so 
that 

I^ujIo (24) 
as done, for example, by iFischer. Henning fc York3 



il994^ . Rigorously speaking, this is true for the first scatter- 
ing only, i.e. when the incoming light is unpolarized. In fact, 
after the first scattering episode, the light acquires some po- 
larization degree and the fraction absorbed by the next event 
does depend on the direction of the electric field with respect 
to the scattering plane, through the following relation 



7 = tj [lo—pi 



1-V 
1 + m' 



(ciQo -I- sif/o)] 



(25) 



which is, in general, different from Ea. l24l Depending on 
the scattering circumstances, this implies that the effective 
albedo can be higher or lower than uj, making this inconsis- 
tent with Eq. 1241 which I have de facto used in the lumi- 
nosity calculations. However, the simulations show that this 
inconsistency applies to photon packets individually, while 
the global effect tends to cancel and the average intensity 
ratio ///o coincides with ui, which can be therefore regarded 
as the average dust albedo (or the albedo for unpolarized 




p/ct 

Figure 13. Radial polarization profiles for a homogeneous sphere 
(ij=50 lyr) with forward (g=0.6, upper) and isotropic (g=0.0, 
lower) scattering, obtained for t=10 yrs. 



light). In my calculations I have used Eq. 1251 even if the 
results are practically identical to those one obtains from 
Eq.m 

The normalized Stokes components after each scatter- 
ing are used as input for the next one, until the packet es- 
capes the dust system at a given projected position on the 
sky plane. 

Since each photon packet is completely independent 
from the others (i.e. there is no phase relation between 
them), the resulting Stokes parameters are additive (CH50). 
This means that the polarization state at any given pro- 
jected position can be obtained simply adding the Stokes 
vectors of all photon packets escaping from that position 
along a given direction. 

In general, multiple scattering tends to act as a depo- 
lari zing mechani sm (see for example ICode fc Whitnevlll995l 
and lLerovlll99jj) . In fact, successive scatterings occur with 
very different and stochastic angular parameters and it is 
very improbable that the polarization gained in the first 
scattering is increased. As a consequence, the maximum po- 
larization ring is expected to become less and less marked, 
making its detection more and more difficult. 

The results for a homogeneous sphere with P=50 lyr, 
which is intended to mimic the case of SN immersed in a 
homogeneous medium, are presented in Fig. ll3l for different 
values for the density in the case of isotropic and forward 
scattering (g=Q.6). The agreement with the analytical so- 
lution is very good for low optical depths. As anticipated, 
for increasing density values the polarization degree tends 
to decrease and the effect tends to be more pronounced for 
forward scattering. 

It must be noticed that in all cases where the dust shows 
an axial symmetry around the line of sight, the net linear 
polarization one would measure in an unresolved LE image 
is null. Conversely, a net polarization degree significantly 
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different from zero would signal the presence of asymmetry 
in the dust distribution, similarl y to the case o f an asymmet- 
ric SN envelope (see for examDle lHoflichl|l99lh . To illustrate 
this effect, I have calculated the net polarization one would 
observe from an unresolved uniform sphere placed at a dis- 
tance D=i?=50 lyr from the SN and seen with a viewing an- 
gle a=0, 7r/4 and 7r/2, which turn out to be <0.1%, 0.8% and 
6.0% respectively, for the usual values of relevant parameters 
{rio — l cm~^, r£j=0.03, (?=0.6, lli—Q.Q). The observed polar- 
ization angle, defined as the direction of the electric field 
vibration, will then give an indication about the direction 
of maximum asymmetry, provided that linear polarization is 
significant with respect to the typical measurement errors. 

Finally, after the second scattering, when the photon 
packets have acquired some degree of linear polarization, 
they may gain also some circular polarization. However, the 
simulations show that it is always Pc <5xl0~*, which is 
within t he noise of our ca l culati ons. Similar results were 
found bv lCode fc WhitnevI (ll995D . For this reason circular 
polarization can be safely neglected (e.g. V—Q). 



12 DISCUSSION 

As I have shown, a LE with an optical luminosity ~10 mag 
fainter than the SN at maximum can be produced by a la 
in the presence of dust with rather normal densities (1-10 
cm~'^) and within 50-100 lyr (see Secj^J. The resulting inte- 
grated light curves have typical slopes of 0.1-0.3 mag yr~^. 
This becomes much higher if the dust is very close to the 
SN, as it is the case, for instance, of a stellar wind. 

For a fixed distance between the SN and a given scat- 
tering element, the contribution to the total flux is not con- 
stant and it depends on its z coordinate (see Fig. , due to 
the presence of forward scattering, which makes the process 
more efficient for foreground dust. Moreover, at any given 
time, the echo luminosity depends on the portion of the 
dusty region intersected by the iso-delay surface. For exam- 
ple, a spherical cloud placed in front of the SN at a distance 
D along the z axis is much more efficient than an identical 
cloud at the same distance D but with its center located on 
the SN plane. 

Of course, the two different geometries would produce 
very dissimilar echo resolved images, but in the lack of LE 
resolved images, it is very difficult to decode the dust ge- 
ometry from the integrated light curve alone. In fact, unless 
the dust is confined within a region with a characteristic di- 
mension 7? <ct, even very different geometries produce quite 
similar light curves (see Secl^J. 

The conclusion is that from the observed integrated 
light curve it is not possible to deduce the dust distribution, 
since the light curve shape is almost geometry independent 
and the echo luminosity can be tuned simply changing the 
dust density and its distance. Therefore, one is left with the 
question whether there are lots of dust far from the SN or a 
small amount of dust close to the SN. 

One may argue that an increase in the dust optical 
depth would also produce observable reddening in the SN 
itself. This is true only under some particular conditions. 
For example, it could well be that there is a distant and op- 
tically thick cloud in front of the SN responsible for a signif- 
icant reddening (but irrelevant for the light echo) while the 



scattering into the line of sight is produced by off-centered 
dust. A good example of such a case is given by SN 1989B, 
which was highly reddened but for which there is no clear LE 
detection (.Milne fc Wel ls 2002). As I will show in a forth- 
coming paper, where I will apply the results presented here 
to the known cases, SN 1998es had a similar behaviour. Due 
to the way the LE forms (Sec.|5J, one can have in fact sev- 
eral possibilities and the LE detection can be accompanied 
by SN extinction or not, in the same way that a heavy SN 
reddening is not necessarily related to the actual onset of a 
LE. 

As a matter of fact, LE properties depend on the global 
dust geometry, while the SN reddening is linked to the dust 
lying on the line of sight only. This is why Td and t^j / have, 
in general, two different meanings, which tend to coincide 
only under special conditions. Therefore, unless one has LE 
resolved images and an estimate of the distance d to the SN, 
the integrated light curve alone is not enough to answer our 
original question, i.e. whether the SN is close or far from a 
dusty region. 

Fortunately, as we have seen, the inclusion of multi- 
ple scattering shows that the LE colour traces pretty well 
the effective dust optical depth (Sec.|7J. So, if an observed 
colour including a reddening sensitive passband (i.e. blue) 
is available, one can draw some conclusions about the dust 
density and the importance of multiple scattering. In fact, 
for low dust optical depths {T^ff 5;0.3), the resulting {B — V) 
colour is similar to the intrinsic colour of the SN at maxi- 
mum light and it is practically independent from the geome- 
try (Sec. l9.ll . As a consequence, any redder observed colour 
would signal the presence of an higher optical depth. In turn, 
this would give an additional constraint when one tries to 
reproduce the observed echo luminosity, thus reducing the 
geometry /density degeneracy which, however, can be elimi- 
nated only when high resolution imaging is also available. 

Under these favorable conditions one can hope to derive 
a 3D mapping of the SN surroundings, given the univocal 
relation between the distance r of the scattering elements 
from the SN and their impact parameter p (see Eq.|SJl. 

It is true that, as MC simulations show, this relation 
tends to break down as the dust optical depth grows, since 
the phenomenon tends to acquire a non-local nature. The 
exact details depend on the dust geometry, but in general 
this takes place at pretty high values of Tg//. In this respect, 
the onset of multiple scattering is much more important for 
the resulting LE luminosity, which does not increase mono- 
tonically as the dust density grows (Sec.[7|l, a fact which is 
foreseen also in the single SSA approximation (Sec.|SJ. Since 
the effect is more pronounced at bluer wavelengths, this pro- 
duces the colour effect I have just discussed. Speaking about 
this, it is worth noting that it is generally assumed that the 
LE has a colour similar to that of the SN at maximum light 
(corrected for reddening if there is any). As I have shown 
(Secl^J, this is surely true when multiple scattering is negli- 
gible but, if the SN is surrounded by a dense dusty medium, 
both the SN and the echo are expected to be reddened, pos- 
sibly by different amounts, since in general Td ^ Te//, so that 
one can have both an unreddened SN coupled to a red echo 
or a re ddened SN associate d with a blu e echo. Both SNe 
1991T JSchmidt et alJll994) and 1998bu iCappellaro et alJ 
l200lf) are examples of the latter instance. In fact, the inte- 
grated [B — V) colour of the LE was about —0.1, while both 
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SNe certainly suffered from significant extinction. The fact 
that the LEs are definitely blue seems to suggest that the SN 
reddening is generated by a fro nt cloud who se characteristic 
dimension R is smaller than ct{2ct + D) (where D is the 
distance of the cloud from the SN itself) , so that it does not 
contribute to the LE for ct > -D + y/D^ + B?, while the 
LE is generated by a less optically thick cloud. I will come 
back to this point in the forthcoming paper. 

Besides the optical depth, which is directly related to 
Cext and plays a relevant role, dust albedo uo and forward 
scattering degree g are also important. For this reason, the 
exact results for LE spectra and broad band colours depend 
on the dust mixture one adopts. As a first approximation, 
in my simulation s I have used t he canonical Ry =3.1 Milky 
Way composition iDraind i2003t) . A possible further develop- 
ment of this work is the study of different dust composition 
effects, similarly to what has been done by Suacrinan ( 2003). 

Of course, the final results are expected to depend also 
on the SN light curve one adopts as input. All the simula- 
tions I have presented in this paper made use of a template 
la input light curve, derived from two normal events (e.g. 
SNe 1992A and 1994D). As I have shown (see Sec.EJ, this 
can be very well approximated by a flash with a duration 
Afsjv, a parameter which normally changes with wavelength 
(see also Table |2J . One may wonder what happens if this 
parameter becomes larger, as it is the case for slow declin- 
ing SNe like 1991T yrs). The answer is easy 
and it is given by Eq. |5] which was derived under simplify- 
ing assumptions: integrated LE luminosity is proportional to 
Afsjv, and therefore slow declining events (e.g. 1991T-like) 
are expected to produce brighter echoes than fast decliners. 
The effect is mild, though, and the total range is of the order 
of 0.5 mag. 

So far I have assumed that the distance d to the SN is 
known, since the main purpose of our inquiry is to derive 
the dust geometry. But the paradigm can be in principle 
reversed, in the sense that the time evolution of a resolved 
LE can be used to infer the distance. This was attemp ted, 
for example, bv lBond etsl\ ^2003^ and lTvlendal J20031 for 
the case of V838 Mon, which has shown an impressive LE. 
This method, in principle promising because it gives a direct 
geometric distance measurement, has nevertheless a weak 
point. In fact, it requires the geometry to be known in order 
to properly use the appar ent ech o expansion to derive the 
distance (the discussion in lTvlenda 2003'). Different geome- 
tries can be of course disentangled looking at the expansion 
evolution, but this requires quite a long time baseline. 

It must be noticed that this technique is applicable 
when the echo edges can be clearly detected and measured 
at different epochs, as in the case of V838 Mon, whose dis- 
tance (6-8 kpc) makes this feasible, while for a SN at a typi- 
cal distance of 15 Mpc this becomes practically not a viable 
method. 

Another application of LE s for direct dis tance determi- 
nation is the one suggested bv lSpark^ il994ft . which makes 
use of the strong polarization ring generated by the dust 
placed at a distance ct from the SN on the z=Q plane (see 
Sec. IIH . Provided that there is dust on the plane of the 
SN, this allows one to get distance estimates directly from 
polarimetric imaging. Since the maximum polarization ring 
obviously expands at the speed of light and not at the su- 
perluminal speed expected for the LE edges, this method is 



practically applicable to old SNe only {t >50 yrs), always 
provided that HST polarimetric data is available. 

There are two problems that could hinder this tech- 
nique. The first is that what counts for the polarization de- 
tection is the polarized fiux which, if any, is smaller than 
the total fiux. As I have shown, for a uniform dust distri- 
bution (which I have mimicked with a sphere centered on 
the SN) the total flux decreases at a rate of several tenths 
of a magnitude per year during the first years, just because 
of the efficiency loss due to forward scattering. This adds 
up to something of the order of several magnitudes in 50 
years, making the detection of the polarized signal very dif- 
ficult, even if its polarization degree remains constant. The 
effect is going to be worsened by the probable presence of a 
bright galaxy background, which would certainly introduce 
additional noise in the polarimetric measurements. Secondly, 
the presence of multiple scattering tends to depolarize the 
scattered light, again decreasing the chances of detection. 
Besides this, I should as well note that a poor imaging res- 
olution tends to reduce t he polarizati on one would actually 
measure (see for example |Ler^^9^ . 

As I have mentioned, the full applicability of this 
method relies on the presence of dust on the SN plane. If this 
is not the case, for example because the region around the 
central source is empty within a radius larger than ct, one 
can obtain only a lower limit for the dis tance. An exampl e 
of such a situation is given by V838 Mon jBond et al .120031) . 
In general, even though the method is certainly promising 
and worthwhile of being pursued, I am persuaded that the 
use of LEs to study the SN environment is more rewarding 
than distance determination. After all, deriving the distance 
to the SN with a factor 2 precision is nowadays not a great 
progress, while the same accuracy on the estimate of the 
distance of a la from a molecular cloud would be of high 
interest, especially if this turns out to be small. 

I like to conclude this discussion with the following con- 
sideration. I have mentioned several times that in pretty 
normal conditions, a la is expected to produce a LE about 
10 magnitudes fainter than its luminosity at maximum. This 
is actu ally what happened f or the two kn own cases of SNe 
1991T JSchmidt et alJll994) and 1998bu iCappellaro et alJ 
I2OOII) which, at around 2 years after the explosion, clearly 
deviated from the usual radioactive decay. As a matter of 
fact, as a resu lt of a sy stemati c LE search inc luding 64 
historical SNe, iBoffi. Sparks fc Macchettd Jl999i) have re- 
ported 16 possible candidates, only one of which is a gen - 
uine la, i.e. SN 1989B (but see also iMilne fc Wellj l20o3l . 
Therefore, one may first inquire why only two events have 
been detected and immediately conclude that this is sim- 
ply because in the vast majority of the cases there is not 
enough dust around la's, not an unexpected conclusion for 
supposedly long-lived and small mass progenitors. 

At least for resolved HST image s, th is point has been 
thoroughly addressed by ISueermanl ||200A in his extensive 
work on the observability of LEs, and his conclusion is that 
many echoes were lost within the surface brightness fluc- 
tuations of background unresolved starlight. In the case of 
ground based observations, I must notice that there are only 
a few la's observed at more than one year past maximum. 
In fact, this has always been a problem, both due to their 
faintness and to the presence of the host galaxy background. 
Therefore, it is difficult to give a final answer on the basis of 
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the current list of LE detections and we will probably have 
to wait a bit more in order to have a statistically significant 
sample. This is even more true if only over-luminous SNe 
are associated with dusty regions. 

What I can exclude here, at least on the basis of my 
calculations, is that very dense (no >30 cm~^), very close 
(r <50 lyr) and very extended {R >30-50 lyr) clouds are 
not present in SNe of type la, for they would produce well 
observable LEs, which is seemingly not the case. Also, from 
the results shown in Sec.0]for the double exponential galac- 
tic disk, it seems that Type la SNe tend to explode far from 
the host galaxy disk, at vertical distances from the galactic 
plane that are significantly larger than the dust height scale 
Zd, otherwise they would produce, again, detectable LEs. 



on the basis of observed LE colors, which can give approxi- 
mate values for re//; 

• The LE spectrum is expected to change very slowly 
with time. It does not resemble any of the SN spectra, being 
a mixture of all SN spectra between t=G and the time of 
observation; 

• Being the product of scattering, a LE is expected to 
be highly polarized, with maximum polarization (~50%) at 
geometrical radius ct. Multiple scattering tends to act as a 
depolarizer, making the maximum polarization ring less and 
less pronounced; 

• In general the polarized flux, which is what counts for 
the polarization detection, decreases with time. This hinders 
the application of polarimetric imaging for direct distance 
determination. 



13 CONCLUSIONS 

In this paper I have discussed the phenomenology of unre- 
solved LEs, with particular attention to the effects of mul- 
tiple scattering on integrated luminosity and broad band 
colours. Even if the treatment is absolutely general and it 
can be applied to LEs produced by any transient source, 
I have devoted my analysis to the SN la case, due to its 
relevance in the possible relation between the observed SN 
properties and the explosion environment. The main results 
I have obtained can be summarized as follows: 

• A type la SN is expected to produce a LE ~10 mag 
fainter than the SN at maximum when it is immersed in a 
medium with dust densities of the order of 1 cm~^; 

• If there is dust on the line of sight, the LE takes place 
immediately after the SN outburst, but it becomes visible 
only when the SN has faded away. This typically happens 
at 2-3 yrs past explosion; 

• In those cases one should rather indicate the phe- 
nomenon as a light reverb to make the parallel with acoustic 
physics more pertinent; 

• In general, the dust optical depth as seen from the SN 
along the line of sight, r^, is not sufficient to describe the 
integrated LE properties. A better description is given by 
Te//, which is a sort of weighted LE optical depth; 

• The integrated LE light curve tends to decline with 
time, both due to geometrical effects and forward scattering. 
Higher forward scattering tends to produce faster declines; 

• For typical values of the forward scattering degree 
(3=0.6) the decline rate during the first few years is of the 
order of 0.1—0.3 mag yrs^^; 

• It is not possible, on the basis of integrated LE curves 
alone, to deduce the dust geometry. This can be inferred 
only if LE resolved images are available; 

• For low values of r^// (<0.3) the LE phenomenon can 
be safely described with single scattering approximation.; 

• For higher optical depths, auto-absorption becomes rel- 
evant and attenuation (re// <1) or, better, multiple scatter- 
ing must be included; 

• Multiple scattering has a number of effects, but the 
most important is the LE reddening, which does not strongly 
depend on the dust geometry; 

• At low dust optical depths (re// <0.3), the LE {B-V) 
colour is similar to the intrinsic colour of the SN at maxi- 
mum. The importance of multiple scattering can be judged 
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